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ABSTRACT 

We present the combination of dynamical accretion model based on 3D GRMHD simulations and 
general relativistic (GR) polarized radiative transfer. We write down the formalism of and perform the 
GR ray-tracing of cyclo-synchrotron radiation through the model of accretion flow in Sagittarius A*. 
GR polarimetric imaging is presented as well as the results for spectrum for a probable set of spins 
and orientations. Precise fitting formulae for Faraday rotation and Faraday conversion coefficients 
are employed for thermal plasma. The axisymmetic flow pattern and the magnetic field geometry 
correspond to averaged 3D GRMHD simulations near the black hole, whereas the analytic model was 
used far from the black hole. The density scaling is found by fitting the sub-mm flux. Spin a = 0.7 
and inclination angle 9 = 0.6 produce the best fit to sub-mm flux and linear polarization fraction. 
Subject headings: 

1. GR POLARIZED RADIATIVE TRANSFER AND DYNAMICS 

We come up with the formalism in several stages. First, we write down the sta.ndard propagation equations of 
Stokes parameters /, Q, t/, V in the uniform thermal plasma |Mclrosc fc McPhedranl[l99l[ l in a locally flat co-moving 
reference frame with synchrotron em issivities/absorptivities from .Melrosel (|1971t) . We take Faraday rotatio n/ conversion 
coeffi cients from lShcherbakovl ()2008f ). as the other published derivation of Faraday conversion coefficients (|Huang et al.l 
l2009t ) is a very crude approximation. Second, we parallel propagate the basis vectors along the null geodesic from the 
observer's plane to account for GR rotation of the basis. Third, with proper gauges on wa ve vector potential we write 
down the covariant equations of polarized radiative transfer. FoUowing lHuang et al.l (120091 ) . we assume that the matrix 
of absorptivities and propagation coefficients generalizes in the polarized transfer analogously to the unpolarized case. 

The dynamical model used in the transfer starts with adiabatic 3D GRMHD simulations of thick accretion flow 
onto the Kerr black hole (BH) with spins a = 0, 0.7, 0.9, 0.98. We average the flow velocity, magnetic field, RMS 
magnetic field, gas density and pressure for the quasi-steady period of the developed accretio n an d separate the electron 
temperature Tg from the proton temperature Tp by applying the heating prescription from lSharma et al.l (j2007l ). The 
dynamical model is smoothly extended to large radii r > 2AM to take into account the Faraday rotation effect at large 
distances from the BH. 
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Fig. 1. — Specific flux Fy in comparison to observations IIYuan. Quataert & Naravanll200^ : IMarronell200"^ (dots) on the left panel. Image 
of specific flux Fi, in logarithmic scale with contrast 8 at = 345 GHz on the right panel. In both calculations spin a = 0.7, inclination 
angle d = 0.6. 
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2. RESULTS 

We perform the GR polarized radiative transfer for each spin for a set of inchnation angles 9 and compare the specific 
fluxes, linear polarization fractions (LP), and circular polarization fractions (CP) to observations. We find, that the 
extreme spins a = 0.9, a = 0.98 do not fit all the observations well. They require lower density n < 2 • 10^ cm^'^ near BH 
to fit the flux at 220 G Hz, but Faraday depo larization fails at these densities leading to high LP at 86 GHz compared 
to the observed valu e ([Macquar t et al.|[2006D . The preferred value of the inclination angle 9 = 0.6 is coincident with 
that in iHuang et al.l ((2009). The spin value a = 0.7 gives the best fit (see Fig. [Ij, though spin a produces good 
fits as well. Imaging produces some unexpected results. As our simulation appears to be substantially sub-Keplerian 
and have significant thermal support, the Innermost Stable Circular Orbit (ISCO) loses its importance and the black 
hole shadow is not always seen. 
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